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1 Introduction 

New phenomena beyond the predictions of the Standard Model (SM) may alter CP 
violation in S-decays. In the decay — )■ J/ip4' the CP- violating phase 0s is the 
weak phase difference between the — B^ mixing amplitude and the b — ?■ ccs decay 
amplitude. SM predicts a small value of 0^ ~ —2/3^ = —0.0368 ± 0.0018 rad [1]. 
Another quantity involved in B^ — B^ mixing is the width difference AF^ = Tl — Th 
of heavy (Bh) and light (Bj^) eigenstates. Physics beyond the SM is not expected to 
affect ATs significantly [2]. Extracting AF^ from data is nevertheless useful as it allows 
theoretical predictions to be tested [2]. The presented analysis uses data collected by 
the ATLAS experiment in pp collisions at a/s = 7 TeV in 2011, corresponding to an 
integrated luminosity of ~ 4.9 fb~^. Measurements of 0^, the average decay width F^ = 
(FL+Fi^)/2 and the value of AF^, were performed with no flavour tagging to distinguish 
between the initial 5° and B^ states. The CP states were separated statistically 
through the time-dependence of the decay and angular correlations amongst the final- 
state particles. 

2 Candidate selection 

ATLAS is a multipurpose detector described in details in [3]. The tracking and muon 
systems are of particular importance in the reconstruction of B mesons. The inner 
tracking detector consists of a silicon pixel detector, a silicon microstrip detector and 
a transition radiation tracker; all surrounded by a superconducting solenoid providing 
a 2 T axial field. The muon spectrometer consists of three superconducting toroids, a 
system of tracking chambers, and detectors for triggering. The triggers used for this 
analysis are based on identification oi a, J/ip ^ At"*"/!" decay, with either a 4 GeV 
transverse momentum (px) threshold for each muon or an asymmetric configuration 
that applies a higher threshold (4 — 10 GeV) to one of the muons. The pairs of 
muon tracks refitted to a common vertex are accepted if the fit results in x^/d.o.f. 
< 10 and the invariant mass falls in the range (2.959 — 3.229) GeV(when both muons 
have \t]\ < 1.05), or (2.852 - 3.332) GeV(both muons within 1.05 < |r/| < 2.5) or 
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(2.913-3.273) GeV otherwise. Candidates for 5° J ii')<p{K+ K-) are sought 
by fitting the four tracks to a common vertex resulting in x^/d.o.f. < 3, while the 
invariant mass of the two muons was fixed to the J/ip mass [4]. In total 131k i?^ 
candidates are collected within a mass range of — )■ K^K~ (1.0085 - 1.0305) GeV 
and 5.15 < m{B^) < 5.65 GeV, figure 1. For each candidate the proper decay time t 
('decay time' henceforth) is determined by the expression: t = L^y M^/{c p^g), where 
is the transverse momentum of the candidate and Mg is the mass of the 5° 
meson (5.3663 GeV) [4]. L^y is the displacement in the transverse plane of the B^ decay 
vertex with respect to the primary vertex (PV) projected onto the direction of Ptb- The 
position of the PV is refitted following the removal of the tracks used to reconstruct the 
Bg candidate. For the selected events the average number of pileup interactions is 5.6, 
necessitating a choice of the best candidate for the P V at which the B^ is produced. The 
variable used is a three-dimensional impact parameter do, calculated as the distance 
between the line extrapolated from the B^ vertex in the direction of the 5° momentum, 
and each PV candidate. The chosen PV is the one with the smallest do- 

3 Analysis and results 

The B^ — J/ipcj) decay involves three angular momentum states of the J/ipcf) system, 
combined into three polarization amplitudes, longitudinal polarization (Aq), and trans- 
verse polarization with the linear polarization vectors of the vector mesons parallel (Ay) 
or perpendicular (A_|_) to each other. The first two states are CP-even, while the last 
state is CP-odd. Another CP-odd state can be produced by a non-resonant K^K~ 
pair or by the decay of the spin-0 /o meson, resulting in another independent ampli- 
tude, the S-wave As- The time evolution of the four decay amplitudes along with six 
interference terms is fitted simultaneously with the angular distributions of the cascade 
decay P° — )■ J/xl:{fi~^fi~)4>{K'^K~). Three independent angles fl = {Ot.iI^t^^t), are 
defined in the transversity basis [5]. The P° J decay probability is a function of 
the following physics parameters of interest, AF^, 0^, F^, the amplitudes Aq, A\\, ^5(0) 
and of the related strong phases defined as 5\\= aig{A\\,Ao), 5± = aig{A±, Aq) and 6s 
= arg(y45, Aq). For untagged analysis all terms involving the mass splitting Am^ in the 
time-dependent amplitudes cancel out [6]. In addition, the time-dependent amplitudes 
depending on 6± are multiplied by a small value of sin 0^, hence the untagged analysis is 
not sensitive to 6±. A Gaussian constraint to the external data, 6± = (2.95 ± 0.39) rad 
[7] is therefore applied. 

An unbinned maximum likelihood fit uses information of the reconstructed P° 
candidate mass m, decay time t, their uncertainties am and at, and the angles Q = 
{9t, i'T, V^t)- In total 26 parameters are determined, including the eight physics param- 
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eters of interest mentioned above, while the other quantities describe the J/ip(f) mass 
distribution, the decay time and the angular distributions of the background. The 
single-event likelihood has the form: 

£ OC ■ /s ■ + /s ■ /b« ■ J^Bf + (1 - /s ■ (1 + /bo)) ■ J'bkgj 

where the index i is used for the variables specific for each single-event, fs is the fraction 
of signal candidates, J^si and J^hkgi are probability density functions (PDF) modelling 
the signal and background. The backgrounds — )> J/ipK* and J/ipK^ir" are 

parametrized separately by J-'^o with /^o being the fraction of this background events. 
The weighting factor Wi accounts for a small decay-time dependency of the acceptance, 
related to a limited resolution in the on-line track reconstruction, details are given in 
[6]. The PDF describing the signal, J-'s, has the form: 

The signal mass density Ps(mj|crmJ is modelled as a Delta function smeared by a 
Gaussian with a per-candidate mass resolution arm ■ Similarly, each of ten terms of the 
signal time and angular dependence, Ps{^i,ti\atJ, is convoluted with a Gaussian with 
a per-candidate resolution at- The angular sculpting of the detector and kinematic cuts 
on the angular distributions is included in the likelihood function through A{Qi,pT^i). 
This is calculated using a four-dimensional binned acceptance method, applying an 
event-by-event efficiency according to the transversity angles {9t, ipT, <Pt) and the pxs- 
The acceptance was calculated from the B^ — )■ J/ ipcf) MC events. The background PDF 
has the following composition: 

^bkgi = Pb(m,) ■ Pb(^mJ ■ Ph{tiWu) ■ PhiOr) ■ Pbi^r) ■ PhM ■ PhK) ■ Pb(PT,) 

The decay time function Pb(^j|o'tJ is parameterised as a Delta function, two positive 
exponentials and a negative exponential. These functions are smeared with the same 
resolution function as the signal decay time-dependence. The prompt peak models the 
combinatorial background events, which are expected to have reconstructed lifetime dis- 
tributed around zero. The two positive exponentials represent a fraction of longer-lived 
backgrounds with non-prompt J/^p, combined with hadrons from the primary vertex 
or from a. B/D hadron in the same event. The negative exponential takes into account 
events with poor vertex resolution. The shape of the background angular distributions, 
-Pb(^T); Phifr), and Pb(V'r) arising primarily from detector and kinematic sculpting are 
described by the empirical functions with nuisance parameters determined in the fit. 
The correlations between the background angular shapes are neglected, but a system- 
atic error arising from this simplification was evaluated. The background mass model. 
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Ph{tTi) is a linear function. Mis-reconstructed — )■ J/tpK* and — J/ipK^ir' 
(non- resonant) decays, are parametrized separately. The fractions of these components 
are fixed in the likelihood fit to values (6.5 ± 2.4)% and (4.5 ± 2.8)%, calculated using 
MC events. Mass and angles have fixed shapes determined from the MC studies. The 
decay time is described by an exponential smeared with per-candidate Gaussian errors. 
Finally, the terms -Ps,b(c"mj5 -Ps,b(o"tJ and -Ps,b(j'Ti) are introduced to account for differ- 
ences between signal and background per-candidate mass and decay time uncertainties 
and values of transverse momenta, details are given in [6]. 

Systematic uncertainties are assigned by considering effects not accounted for in the 
likelihood fit. The impact of inner detector residual misalignments was estimated using 
events simulated with perfect and distorted geometries. Systematics due to limitations 
of the fit model were determined by 1000 pseudo-experiments generated with variations 
in the signal and background mass model, resolution model, background lifetime and 
background angles models. Systematics due to B^ — )■ J / tpK*^ and B^ — )■ J/ tpKir arise 
from the uncertainties of the PDG decay probabilities, ref. [4]. 

In the absence of initial state flavour tagging the PDF is invariant under the simul- 
taneous transformations: {0s, AF^, S±, S\\, Ss} {—(ps, AF^, vr — S±, —S\\, —6s} leading 
to a fourfold ambiguity. As the constraint on 6± is taken from the LHCb measure- 
ment [7], that quotes only two solutions with a positive (ps and two AF^ values sym- 
metric around zero, two of the four minima fitted in the present non-flavour tagged 
analysis are excluded from the results presented here. Additionally a solution with 
negative AF^ is excluded following the LHCb measurement [8] which determines the 
AFj, to be positive. The measured values, for the single minimum resulting from these 
constraints, are given in Table 1. The second strong phase, 6\\, is fitted very close to 
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Value 


Stat. 


Syst. 


Par. 


Value 


Stat. 


Syst. 


0s(rad) 


0.22 


0.41 


0.10 


\M0)\' 


0.528 


0.006 


0.009 


AF,(ps-i) 


0.053 


0.021 


0.010 


l^||(0)P 


0.220 


0.008 


0.007 


rs(ps"^) 


0.677 


0.007 


0.004 


\M0)\' 


0.02 


0.02 


0.02 



Table 1. Fitted values for the physics parameters along with their statistical and systematic 
uncertainties. 

its symmetry point at tt. Pull studies, based on pseudo-experiments using input values 
determined from the fit to data, return a non-Gaussian pull distribution for this pa- 
rameter. For this reason the result for 6\\ is given in the form of a la confidence interval 
[3.04,3.24] rad. The strong phase of the S'-wave component is fitted relative to 6±, as 
6± — Ss = (0.03±0.13) rad. The fraction of S'-wave KK or /o contamination is measured 
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Figure 1. Mass and decay time fit projections for the candidates. The pull distribution 
at the bottom shows the difference between the data and fit value normalised to the data 
uncertainty. 



to be consistent with zero, at |y45(0)p = 0.02 ± 0.02. The two-dimensional likeUhood 
contours in the (ps — AFg plane for the 68%, 90% and 95% confidence intervals are 
produced using a profile likelihood method and are shown in figure 2. The systematic 
errors are not included in figure 2 but as seen from table 1 they are small compared 
to the statistical errors. The ATLAS measured parameters of the — )■ J/ip4> decay, 
using 4.9 fb~^ of integrated luminosity collected in 2011, are consistent with the world 
average values and with theoretical expectations, in particular (pg is within la of the 
expected value in the Standard Model. 
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5 1 constrained to 2.95 ± 0.39 rad — 68% C.L. 



- 90% C.L. 
■ 95% C.L. 
■- Standard Model 

Ar, = 2|r,jcos((^j 




Figure 2. Likelihood contours in the (ps — AF^ plane, statistical errors only. The green 
band is the theoretical prediction of mixing-induced CP violation. 
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